NASA Technical Memorandum 101417 
AIAA-89-0437 


Numerical Analysis of Flow Through 
Oscillating Cascade Sections 


{NASA-'iM— 101^*17) SiUMEBICAL 
THROUGH CSClLLATiEG CASCADE 

p 


A E AC YS IS OF FLOW 
SECIICKS (MASA) 
CSCL 01 A 


G3/02 


N89- 14220 


Unclas 

0179840 


Dennis L. Huff 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

27th Aerospace Sciences Meeting 

sponsored by the American Institute of Aeronautics and Astronautics 
Reno, Nevada, January 9-12, 1989 





NUMERICAL ANALYSIS OF FLOW THROUGH OSCILLATING CASCADE SECTIONS 


Dennis L. Huff 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Abstract 

A compressible, Euler or full Navicr-Stokcs, finite 
difference code has been developed for modeling 
inviscid and viscous flow through oscillating cascades. 
The code uses a deforming grid technique to capture 
the motion of the airfoils and can model oscillating 
cascades with both zero and non-zero inter-blade 
phase angles. Two-dimensional, unsteady character- 
istic boundary conditions are applied at the inlet for 
viscous solutions and to both the inlet and exit for 
inviscid solutions to minimize reflections from these 
regions. 

Results are presented for two cascade geometries 
for comparisons with experimental data and a flat 
plate cascade for comparisons with small- 
perturbation theory. Overall, there is good agree- 
ment for the two non-zero inter-blade phase angle 
cases presented and poor agreement for the zero 
inter-blade phase angle cases. Studies on reflecting 
versus non-reflecting inlet and exit boundary condi- 
tions show that the treatment of the boundary can 
have a significant effect on the first harmonic, 
unsteady pressure distributions for the cases with 
zero inter-blade phase angles. Using first or second- 
order temporal accuracy in the numerical algorithm 
did not make a significant difference in the unsteady 
pressure distributions for the present solutions. 


Nomenclature 

C sonic velocity 

CP pressure coefficient, 





c blade chord length 

e total energy of the 

fluid per unit volume 

g/c gap-to-chord ratio 


i 

incidence angle 

lm{ } 

imaginary part of { } 

J 

Jacobian of 
transformation 

k 

reduced frequency based 
on semichord, j-r 

M 

Mach number 

P 

total pressure 

P 

static pressure 

Rc 

Reynolds number based 
on chord 

Re{ } 

real part of { } 

s 

arclength of a grid 

line in the 
redirection 

t 

time normalized by ^ 

u,v 

Cartesian velocities 
normalized by C x 

V 

total velocity 

w 

weighting function for 
grid deformation 


Cartesian coordinate 
system for cascade 

* ' . y ' 

Cartesian coordinate 
system for airfoil 

a 

amplitude of pitching 
motion 

P 

flow angle 

AC P 

pressure difference 

p. p. 

coefficient, , 
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ACP' pressure difference 

/j , ' p . 

coefficient, — — 

p , V \ i a| 

Y stagger angle 

*1 normal direction of 

transformed coordinate 
system 

^ chordwise direction of 

transformed coordinate 
system 

4> phase lead angle 

relative to the 
complex angular 
displacement 

p fluid density 

t time variable 

c ' thickness-to-chord 

ratio 

v inter-blade phase angle 

co airfoil oscillation 

frequency 

1 conditions at the inlet 

2 conditions at the exit 

+ upper and lower 

surfaces on airfoil, 
respectively 


Introduction 

The design of turbomachinery blades requires the 
prevention of flutter for all operating conditions. 
However, flow field predictions used for aeroelastic 
analysis are not well developed for all flow regimes. 
Many methods used for flutter analysis depend on the 
successful predictions of blade loading and blade 
motion. Ideally, this is an interactive process where a 
structural analysis determines the blade motion from 
the blade loading and the loading is determine by the 
flow analysis from the blade motion. Unfortunately, 
the flow analysis technology lags the structural analy- 
sis technology in that accurate predictions of the flow 
field for transonic flow through oscillating cascades 
are not readily available. Many of the existing flow 
analysis methods are numerical and require large 


amounts of computer time. However, advances in 
computer technology have given the aerodynamicist 
hope for numerical methods being used as a tool in 
the design process. 

In theory, the solution of the Navier-Stokes 
equations should model all flow phenomena asso- 
ciated with transonic flow through oscillating cas- 
cades. Since exact solutions arc not available, 
computational techniques are used and introduce 
approximations to the solutions. These approxima- 
tions commonly use eddy viscosity methods for turbu- 
lence modeling, numerical dissipation terms for 
stability and coarse grids for practical computing. 
Further complications arise when unsteady flow exists 
and proper treatment of the boundary conditions are 
essential. The message here is that numerical solu- 
tions of the more complex flow equations, like the 
Navier-Stokcs equations, do not automatically give 
accurate results. Much research is needed in this 
area to develop proper modeling of the physical flow 
field. 

The present research focuses on numerical 
solutions of the Euler and Navier-Stokes equations to 
model two-dimensional flow through oscillating cas- 
cades. The model prescribes harmonic pitching 
motions for the blade sections for both zero and 
non-zero inter-blade phase angles. The deforming 
grid technique introduced in reference 1 is utilized 
for convenient specification of the periodic boundary 
conditions. Several sample predictions are compared 
to experimental data. Also, an investigation of 
reflecting versus non-reflecting boundary conditions 
is presented to determine their effect on the unsteady 
blade loading. 

Governing Equations 

A major portion of the present code is based on 
the unsteady, viscous solver developed by Sankar and 
Tang (ref. 2) for flow past isolated airfoils. This code 
solves the two-dimensional, unsteady, Rcynolds- 
averaged, compressible Navier-Stokes equations on a 
body-fitted moving coordinate system in strong 
conservation form using an ADI procedure. These 
equations can be written as; 

f /< + V + 5 n (1) 
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where 


q = J~'{ p.pit.pi'.e} 

and pis the fluid density; u and v are the Cartesian 
components of the fluid velocity; c is the total energy 
per unit volume. The body-fitted (?, ■ *1 • t ) coordinate 
system is related to the Cartesian coordinates using 
the following transformation: 

$-$(*. y.O 

T| = ii(x,y.t) 

T * t 


The Jacobian of the transformation is given by: 


J = ^xll y -bx f .y 


1 

V|y, 1 -v 11 y l 


CO 


Grid 

An unique feature ol the present code is the 
treatment of the grid for oscillating cascades. A 
method for deforming the grid was developed in ref- 
erence 1 for zero inter-blade phase angles. The 
present study is an extension of this method to model 
non-zero inter-blade phase angles. ^ 

The code uses a C-grid generated from Sorenson’s 
(ref. 7) GRAPE code, which was modified by Chima 
(ref. 8) for improved modeling in turbomachincry 
problems. One C-grid is generated for each blade in 
the cascade. The outer boundary of the C-grid is 
defined by the user in the GRAPE code. A deform- 
ing grid technique (ref. 1) is used to locate the posi- 
tion of the grid as a function of time. The inner 
boundary moves with the prescribed blade motion, 
while the outer boundary remains fixed in space. The 
grid lines connecting the inner and outer boundaries 
are allowed to deform. The amount of deformation 
is a function of the distance away from the surface of 
the airfoil. A weighting function, w, is defined as: 


and the metrics of the transformation arc given by: 




Standard central differences were used to compute 
x t . y t . -X,,, y„ which were then used to calculate the 

metrics. . . 

The P and 0 terms in equation (1) are the inviscid 

terms in the \ and redirections, respectively. The 
viscous terms, k and 5, are treated explicitly and can 
be omitted to give solutions for the Euler equations. 
The Beam- Warming, block ADI algorithm is used to 
solve the governing equations. Artificial dissipation is 
added to help the stability. The solution is second- 
order accurate in space and first or second-order 
accurate in time. The Baldwin-Lomax, two-layer 
algebraic model (ref. 3) is used to evaluate the eddy 
viscosity. The Johnson-King (ref. 4) and Gorski (ref. 
5) Jt - c turbulence models have been coded by Wu 
(ref. 6), but no comparisons arc included in this 
study. Further information about the algorithm for 
isolated airfoils can be found in reference 2. 


q) - 


s(5.n) 

S ( . Umax) 


( 6 ) 


where s = the arclength of a grid line from the airfoil 
surface (n= 1) to some grid point along f, = constant, 
and it,,,,, = the outer boundary grid line. The grid 
deformation is defined as: 

A x (< = w„(Ax„) 

Ay,/ “ ( Ay ,y ) ( 7 


where Ax,,and Ay u arc the spatial differences 
between successive time steps if the entire grid was 
moved as a rigid body. From equations (6) and (7), 
we sec that nodes at the inner boundary (s=0) gives 
Wj) - | ? which means the airfoil surface follows the 
rigid body motion of the blade. Conversely, the outer 
boundary nodes give 0 and the node positions 
remain fixed at the initial specified locations. The 
interior nodes shear in space relative to the initial 
grid as w„ varies between 0 and 1. The node veloci- 
ties can be easily found by dividing the grid deforma- 
tion by the lime step value. 

Multiple blade computations arc made possible by 
stacking the C-grids for each blade and passing infor- 
mation between the upper and lower boundaries. 
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Each C-grid is expanded by one grid line in (he 
redirection at the outer boundary to provide ghost 
points for the interior grid of the adjacent blades. 
This allows the periodic boundary condition to be 
treated implicitly. The case for zero inter-blade 
phase angle is the simplest for grid generation. Peri- 
odic boundary conditions are applied across the 
upper and lower boundaries thereby requiring a grid 
for only one blade. However, non-zero inter-blade 
phase angles require multiple blade computations for 
exact treatment of the periodic boundary conditions. 
For example; an inter-blade phase angle (a) of ±180 
degrees requires two blades, o = ±90 degrees 
requires four blades, etc. This significantly increases 
the computational time for non-zero inter-blade 
phase angles, but provides an exact boundary condi- 
tion. This method may be useful as a tool for evalu- 
ating simplified periodic boundary conditions that do 
not use as many blades for non-zero inter-blade 
phase angles. 

The deforming grid technique is demonstrated in 
Figure 1 for a multiple blade computation. The 
figure shows two extremely simplified grids superim- 
posed at two different times. Multiple grids arc 
generated for a NACA 16-004 cascade with twenty 
degrees of stagger. This geometry is representative 
of a geometry that might be used in the advanced 
turboprop. The inter-blade phase angle is 90 degrees 
and requires four blades to specify the periodic 
boundary conditions. (Only three blades are shown 
in this figure for clarity.) The amount of oscillation is 
exa 8& cr ated to illustrate the deforming grid tech- 
nique. Notice how the outer boundary of the grid 
around one blade remains fixed in space, while the 
inner boundary follows the motion of the airfoil. An 
actual calculation for an oscillating cascade typically 
has a pitching amplitude of only a few degrees and 
does not distort the grid as much as shown in Figure 
1. Larger pitching amplitudes are possible, however, 
by regriding at different time periods and using the ’ 
deforming grid technique for various time intervals of 
the pitching motion. 


The proper treatment of boundary conditions is 
important for any numerical computation and espe- 
cially for unsteady flows. The present solution solves 
for the flow around each blade independently and 
uses periodic boundary conditions along the upper 
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and lower boundaries to model the cascade effects. 
Ghost points are assigned at the first interior grid line 
( ] 1 ~ fimax “ I) and arc used implicitly by the adjacent 
grid from the next blade. Although it is tempting to 
use the most current flow information as it becomes 
available from the integration scheme, it is important 
to only use flow information from the same time 
period across the periodic boundaries. This elimi- 
nates time inaccuracy due to the direction of time 
marching from blade to blade. The metric data is 
also forced to be continuous along the periodic 
boundaries. This procedure essentially makes the 
periodic boundaries invisible to the flow solution and 
means the only specification of boundary conditions 
are at the inlet and exit planes, blade surfaces and the 
slits aft of the airfoils. 

The inlet conditions are assumed to be uniform 
with a reference density, velocity and inlet flow angle. 
In previous solutions (ref. 1), the static pressure was 
also specified from the steady-state solution of a sta- 
tionary cascade. This provided enough information 
to determine the unsteady solution for an oscillating 
cascade, but also introduces a reflective boundary 
condition. This problem has been addressed by 
transforming the fluid equations into characteristic 
form at the inlet boundary and prescribing the flow 
information as incoming or outgoing waves. This 
approach was formulated by Hedstrom (ref. 9) for 
one-dimensional gas dynamic equations, and 
extended to two-dimensional problems by Thompson 
(ref. 10). While the two-dimensional method does 
not give a perfectly transmitting boundary condition, 
it does offer a method for minimizing reflections by 
considering waves oblique to the grid boundary. A 
brief explanation of this method is given below. 

The two-dimensional fluid equations from equa- 
tion (1) can be generalized in the following form if we 
neglect viscous terms: 


*0 dP dC 

dt + ^ + <jri + C< + C '’ = 0 


Define the following relations: 


dO du dP du 

dt dt d% V 


dC „ dU 

j^-R-.A = P- Q.B-P 'R, 


( 8 ) 


4 



where 


c % -r'c x .c n =p 'c n (9) 

Substituting into equation (8) and multipying by P 
gives: 


bU 


■V 


\i k m 


k "' k br] 
0 


for outgoing waves 
for incoming waves 


bU_ 

bt 


dll J U r- 


( 10 ) 


The 71 and B matrices can be defined as: 

/t = s' 1 AS,fl=r l Hr (H) 


where A and H are diagonal matrices which define 
the characteristic velocities. Hence, the characteristic 
form of equation (8) can be written as: 


Here, m* are the left eigenvectors of B, u* are the 
eigenvalues of B, and — is evaluated using one-sided 
differencing. The problem now becomes solving for 
— from equation (16) and computing '77 from equa- 

tion (13): 


dO dP ■ du 

— + — + C. - P — 

di 5 a/ n 


(17) 


d JL 

Jt 


bU , &U 
'AS-^ + T'ur — 
a^ aq 


+ c ? + c n =o 


(17) 


A portion of the q ma x grid l' nc ‘ s defined to be the 
inlet to the cascade. This requires equation (8) to be 
written in characteristic form only in the redirection: 

dO dP 

a/ d% * 

+ p ( t ’' HT ^ +c n) = o (13) 


Consider the terms in the parentheses containing the 
information in the indirection. Define: 


' , dU „ 

r Hr — + c, 


-if <“> 


Multiply equation (14) by T: 


T ir +H7 r + 7C ’ =0 

a< n aq 


(K>) 


The boundary conditions at the inlet boundary can be 
written in a form similar to equation (15): 


dU_ 

l ‘at„ 


+ n t + m t c„ = 0 . 


(16) 


The boundary conditions applied at the inlet for sub- 
sonic inflow specify density, velocity and flow angle 
and solve equation (17) to determine the ener©'. 

For viscous flows, the exit boundary conditions 
extrapolate density and velocity from the interior and 
specify static pressure to calculate energy. The char- 
acteristic equations are not valid across the viscous 
wake due to the non-isentropic flow at the exit 
boundary. This introduces reflective boundary condi- 
tions at the exit. It is not clear how to apply a 
non-reflective boundary condition in this region. For 
inviscid flows, the characteristic equations are used by 
specifying static pressure at the exit and solving an 
equation similar to equation (17) for density and the 
two components of velocity. 

Solid wall boundary conditions are applied along 
the airfoil surface and the flow variables are averaged 
across the slit aft of the airfoil. 

and Discussion 

NACA 65-S eries Cascade 

Sample predictions were presented in reference 1 
to validate the use of a deforming grid for an isolated 
airfoil. The present study investigates a NACA 65-se- 
rics cascade with subsonic inflow to help validate the 
solver for both zero and non-zero inter-blade phase 
angles. The predictions presented for this cascade 
are obtained from the solution of the full Navier- 
Stokcs equations and therefore include viscous losses. 
The experimental data was obtained by Carta (ref. 

11) and has been used by other researchers for code 
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validation. This cascade has been called "The First 
Standard Configuration” by the Second International 
Symposium on Aeroelasticity in Turbomachines (ref. 
12). Figure 2 shows the cascade geometry nomencla- 
ture used in the present investigation. 

The cascade geometry consists of NACA 65-series 
thickness distributions on a 10 degree circular-arc 
camber line, 55 degree stagger angle (y), a gap-to- 
chord ratio (g/c) of 0.75 and a thickness-to-chord 
ratio (t ') of 0.06. An arbitrary test case from the 
First Standard Configuration has been chosen for 
comparison with the present solver. The inlet Mach 
number (Af x ) is 0.17 and the airfoils arc pitching 
about the midchord with an amplitude (a) of 2.0 
degrees and a reduced frequency (k) of 0.123 based 
on semi-chord. Two inter-blade phase angles (a) of 0 
and 90 degrees are investigated in this study for com- 
parisons with experimental data. 

The grid used in the solution for this cascade is 
shown in Figure 3. The GRAPE code is used to 
generate a grid about the mean pitching angle. The 
solver is used to determine the mean flow solution. 

An iterative process is necessary by the user to verify 
that the specified exit pressure gives the desired 
velocity specified at the inlet. (Losses through the 
cascade prevent an exact analytical specification of 
these properties.) Four blades arc required for o=90 
degrees, although only two arc shown in Figure 3a. 
Multi-blade solutions for oscillating cascades are 
done by generating grids for each blade and using the 
mean flow solution as an initial condition to the 
unsteady solutions. The solver automatically gener- 
ates these grids before the unsteady solution begins 
by deforming the mean flow grid for one blade 
through one cycle of oscillation and saving the grids 
that occur at multiples of the desired inter-blade 
phase angle. For example, when a=90 degrees, four 
grids need to be assembled with a blade-to-blade 
phase angle of 90 degrees. The grid from the mean 
flow solution is used for the first blade (<jU = 0). The 
grids for the remaining three blades arc found by 
deforming the mean flow grid through one cycle of 
oscillation and storing the grid coordinates when 
ca1 = 90, 180, and 270 degrees. The solver numbers 
the blades and essentially stacks the grids to give one 
global grid containing four blades. The global grid is 
then used as the initial grid for the unsteady solution, 
where it is deformed for the oscillating cascade using 
the method described in the "Grid” section. Since the 
dynamic memory is usually limiting for multi element 


grids, the solid-state storage device on the CRAY- 
XMP at NASA Lewis Research Center is used to 
store all grids and flow information. The C-grid and 
flow variables for one blade are transferred into 
dynamic memory on the CRAY as needed by the 
solver. 

Figure 3b shows more detail of the grid near the 
airfoil surface. The distance of the first grid line off 
the airfoil surface is 0.00005 chord lengths, which is 
appropriate for modeling the boundary layer. The 
grid around one airfoil is 157 x 40 in the ^ and 
M-directions, respectively. The grid is clustered about 
the leading edge and the boundary layer region. 

Also, the grid near the inlet is clustered in the redi- 
rection to resolve waves near the inlet. 

A steady-state solution was done first for the 
NACA 65-series cascade for the mean flow condi- 
tions. This required adjusting the exit pressure until 
the inlet Mach number was near 0.17. As mentioned 
by Verdon (ref. 13), there is some ambiguity concern- 
ing the inlet flow angle. The experimental flow angle 
was P i = 66.0 degrees. Verdon found that P , = 62.2 
degrees gave better agreement with the experimental 
mean flow pressure distributions. A study of mean 
flow pressure distributions for various inlet flow 
angles is presented in Figure 4 using the present 
solver. An inflow angle of 64 degrees was found to 
give the best comparison with the experimental data. 
Each solution required 2000 time steps to reach a 
steady-state and used about 400 seconds of CRAY- 
XMP CPU time. 

Two unsteady solutions for inter-blade phase 
angles of 0 and 90 degrees were done using the mean 
flow solution as an initial condition. Both solutions 
were run with first-order temporal accuracy. The 
surface pressure time histories were recorded and 
found to reach a reasonably periodic solution after 
two cycles of cascade oscillation. A Fourier trans- 
form was done on the second cycle to determine the 
first harmonic magnitude and phase of the pressure 
distribution relative to the airfoil motion. Higher 
harmonics of the surface pressure distributions were 
found to be small compared to the fundamental fre- 
quency. The pressures are normalized by the airfoil 
pitching amplitude and the phase is referenced to the 
airfoil pitching angle starting at the maximum (nose 
up) blade angle. The predictions for o= 90 and 0 
degrees are shown in Figure 5 along with the exper- 
imental data. Both the magnitude and phase are well 
predicted for both the upper and lower surfaces. 
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The phase angles for o=0 degrees differ quantita- 
tively on the upper surface toward the trailing edge, 
although the overall trend is predicted. The reason 
for this is probably due to the pressure magnitudes in 
this region being close to zero. The unsteady calcula- 
tions of phase are more sensitive when the output of 
the real and imaginary parts of pressure from the 
Fourier transform go towards zero. The unsteady 
solutions use 3.26x 1 0" 5 seconds of CPU per time 
step per grid point per blade and required 16053 time 
steps to complete 2.25 cycles of oscillation. 

Flat Plate Cascade 

A cascade of flat plates was used in the present 
analysis for comparisons with the unsteady, small- 
perturbation, subsonic analysis in reference 16. The 
cascade was arranged with v = 53 degrees and 
g/c= 0.767, which corresponds to the conditions used 
in the biconvex airfoil cascade to be presented in the 
next section. All runs for the remaining solutions are 
inviscid so that the non-reflecting boundary condi- 
tions are valid at both the inlet and exit planes to the 
cascade. This also allows comparisons between runs 
using reflecting and non-reflecting boundary condi- 
tions. 

The first set of flow conditions consider M ,-0.65, 
p ,=53.0, k= 0.221, and a=0.10 degrees. This gives a 
zero mean incidence angle on the flat plates and a 
small amplitude of oscillation to allow for compari- 
sons with the small-perturbation, subsonic flat plate 
analysis of Smith (ref. 16). A 199 x 33 grid was used 
around a flat plate with x =0.005. The leading and 
trailing edges were rounded to aid the C-grid genera- 
tion, and therefore this is only an approximate repre- 
sentation of a flat plate. Figure 6 shows the 
comparisons for o= -90, 90, and 0 with the 
small-perturbation theory. (Notice that the pressure 
coefficient is now normalized by p , V \ instead of 
ip , V \ to be consistent with the normalization used 

by Smith in reference 16). Also shown in Figure 6 
are comparisons from the present analysis for reflec- 
tive versus non-reflective boundary conditions. The 
"reflective boundary conditions" means the inlet and 
exit planes specify the information from the mean 
flow solution. The "non-reflective boundary condi- 
tions" use the characteristic method described earlier 
to help minimize reflections from the inlet and exit. 

For o=-90 and o= 90, the agreement is very good 


for all cases and the type of boundary conditions had 
little effect on the results. When o= 0, the agreement 
with small-perturbation theory became worse and the 
type of boundary conditions used at the inlet and exit 
gave different the results within the code. In particu- 
lar, the magnitude distribution dropped to values 
closer to the theory when using the non-reflective 
boundary conditions, but the phase values became 
worse. The reason for this is unknown, although this 
condition is very close to resonance, as calculated by 
Verdon (ref. 13). Acoustic resonance occurs in cas- 
cades when a wave propagates from an airfoil in the 
direction of the stagger line and intersects an 
adjacent blade in the same amount of time associated 
with the inter-blade phase angle. Numerical solu- 
tions are expected to be difficult when operating near 
an acoustic resonance condition. 

Solutions were also done for oscillating flat plates 
with M ,=0.80 and k=0.185 and are compared to the 
small-perturbation theory. The results are presented 
in Figure 7 and show the same trends reported for 
Figure 6. The non-zero inter-blade phase angle solu- 
tions are in good agreement with the small- 
perturbation theory. The near-resonant condition , 
o=0, shows poorer agreement with the theory. In 
general, the non-reflecting boundary conditions show 
good agreement in the magnitudes of the unsteady, 
pressures when compared to the theory, but differ in 
the phase distributions. On the other hand, the cases 
using reflective boundary conditions show worse 
agreement in magnitude predictions and better 
agreement in phase distributions. 

NASA Le wis Cascade 

Another cascade geometry for which experimental 
data (including unsteady surface pressure distrib- 
utions) exists is the NASA Lewis Transonic Oscillat- 
ing Cascade (ref.14 and ref. 15). This cascade 
consists of symmetric biconvex airfoils (t =0.076) 
arranged with y= 53 degrees and g/c = 0.767. Both 
subsonic and transonic test cases have been chosen 
from the experimental data for comparisons with pre- 
dictions from the present solver. Again, all runs for 
the remaining solutions are inviscid so that a 
boundary condition study can be performed. 

The subsonic test cases consider M , = 0.65 and the 
mean p ,=60.0 degrees (i=7 degrees). The grid gen- 
eration is similar to the technique described for the 
NACA 65-series cascade. The distance of the first 
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grid line off the airfoil surface is 0.01 chord lengths. 
The grid around one airfoil is 199 x 22 in the f, and 
redirections, respectively. Unsteady solutions were 
done for an oscillating cascade with the biconvex 
airfoils, ot = 1.2 degrees, k =0.221, and o=-90,90 and 0 
and arc shown in Figures 8,9, and 10, respectively. A 
Fourier transform was done on the third cycle to 
determine the first harmonic unsteady pressure dis- 
tributions. For completeness, both the CP and A CP 
predictions are presented for real and imaginary 
parts, and magnitude and phase components. The 
cases for o=- 90 (Figure 8) and «=90 (Figure 9) are in 
good overall agreement with the experimental data of 
reference 15. When o= 0 (Figure 10), there is a shift 
in the real part (in-phase part) of the unsteady pres- 
sure on both the upper and lower surfaces relative to 
the experimental data. Again, operating near an 
acoustic resonant condition may contribute to these 
discrepencies. It is interesting to note that while the 
upper and lower unsteady surface pressures arc in 
poor agreement with the data, the A CP representa- 
tion of the same predictions look better. This illus- 
trates the importance of presenting separate pressure 
distributions for the upper and lower surfaces when 
possible. 

An investigation of the numerical time accuracy is 
done for the same test cases presented above. The 
code can be run with cither first-order or second- 
order accuracy in time. Figure 11 compares the real 
and imaginary parts of pressure for o=-90, 90, and 0, 
and shows that there are only minor differences 
between the predictions for first and second-order 
temporal accuracy. 

The transonic test cases arc identical to the 
subsonic test cases, except M , = 0.80. The increase in 
the inlet velocity causes a shock to form on the upper 
surface near the leading edge. Figure 12 shows the 
mean flow pressure distributions for two conditions: 
l.) matching the experimental inlet Mach number 
(0.80) and 2.) matching the experimental ratio of the 
exit static pressure to the inlet total pressure (0.7248). 
The desired flow conditions probably lie somewhere 
between these predictions, as shown by the compari- 
son with experimental data. In both cases, the 
pressure distributions are in good agreement with 
each other and the experimental data and therefore 
either condition could be used for the mean flow in 
the unsteady analysis. The unsteady solutions (that 
follow) match the inlet flow velocity from the exper- 
imental data. 


The unsteady pressure distributions for k = 0.185, 
°= -90,90, and 0 degrees arc shown in Figures 13, 14, 
and 15, respectively. The change in the reduced 
frequency maintains the same oscillation frequency 
used in the subsonic cases ( oo=200 Hz). The agree- 
ment with experimental data are very good for the 
o=-90 degrees case (Figure 13). The shock on the 
upper surface is modeled by the analysis, as evident 
by the pressure peaks near the leading edge. Note 
that the experiment did not have pressure transducers 
near this region and is not expected to capture this 
behavior. The agreement with the data becomes 
progressively worse for o= 90 (Figure 14) and o=0 
(Figure 15). The real-part of pressure on the forward 
portion of the upper surface differs from the data for 
ct=90. This was also observed by Verdon (ref. 13) 
using a linearized potential code. Both the real and 
imaginary parts of pressure differ from experiment 
for o=0. Again, this case is near an acoustic reso- 
nance condition. Also, the real and imaginary parts 
of pressure arc significantly smaller than other cases, 
which makes the calculations for phase sensitive to 
their values. 

A study was done to determine the effect of 
boundary conditions on the first harmonic pressure 
distributions. Each of the three cases reported above 
were run with four combinations of inlet and exit 
boundary condition types: 1.) reflecting inlet and exit, 
2.) "non-reflecting" inlet and exit, 3.) reflecting inlet, 
"non-reflecting" exit, and 4.) "non-reflecting" inlet, 
reflecting exit. The results are presented in Figure 16 
and show that the type of boundary conditions can 
make a significant difference on the first harmonic 
unsteady pressure distributions. For clarity, the indi- 
vidual cases are labeled only when there are noticable 
differences in the pressure distributions. The type of 
boundary condition used has the little effect when 
°= -90 degrees. However, when o= 90 degrees, the 
inlet boundary condition has a significant effect, par- 
ticularly near the leading edge shock. The zero inter- 
blade phase angle case, which is near an acoustic 
resonance condition, shows a different solution for 
each type of boundary condition. This is not to say 
that these observations will hold true for other flow 
conditions. Much work has been done by other 
researchers to predict when acoustic waves will prop- 
agate from an oscillating blade row. They have 
shown that different geometries and flow conditions 
will change the propagation characteristics. While a 
full study of wave propagation is beyond the scope of 
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this work, the present study does show that proper 
treatment of the inlet and exit boundary conditions is 
neccessary for modeling the first harmonic unsteady 
pressure distributions in transonic flows. Overall, the 
results from the boundary condition studies show that 
the type of inlet boundary conditions have a greater 
effect on the unsteady pressure distributions than the 
type of exit boundary conditions. This is encouraging 
for validating the viscous solutions that cannot use 
the non-reflective boundary conditions at the exit. 
However, this conclusion is not substantiated without 
a further study on the effects of inlet and exit distance 
from the blade surface. All of the present cases have 
an inlet boundary closer to the blade than the exit 
boundary. Therefore, the downstream-running waves 
have a longer distance for attenuation than the 
upstream-running waves. 

The inviscid runs require about ?. 3b.v 1 0 
seconds of CPU per time step per grid point per 
blade. The number of time steps required for one 
cycle of oscillation is a function of k,Af h and A t. The 
cases with k = 0.185, M , = 0.80 and At = 0.01 require 
2835 seconds of CPU for 6899 iterations, which cor- 
responds to 3.25 cycles of oscillation of four blades. 

It is possible to reduce this run time by performing a 
grid and time step size study. 

Conclusions 

A compressible, Euler or full Navicr-Stokes, finite 
difference code has been developed for modeling 
inviscid and viscous flow through oscillating cascades. 
The code uses a deforming grid technique to capture 
the motion of the airfoils and can model oscillating 
cascades with both zero and non-zero inter-blade 
phase angles. Two-dimensional, unsteady character- 
istic boundary conditions are applied at the inlet for 
viscous solutions and to both the inlet and exit for 
inviscid solutions to minimize wave reflections from 
these regions. Results show that predictions for a low 
speed, NACA 65-series oscillating cascade are in 
good agreement with experimental data. The predic- 
tions for an oscillating cascade of flat plates are in 
good agreement with small-perturbation theory for 
non-zero inter-blade phase angles. The zero degree 
inter-blade phase angle cases, which were near an 
acoustic resonant condition, difter from the theory. 
The predictions for an oscillating cascade of biconvex 
airfoils are in fair agreement with experimental data 
for non-zero inter-blade phase angles. Solutions for 


zero degree inter-blade phase angles only show quali- 
tative agreement with the data. Studies on reflecting 
versus non-reflecting inlet and exit boundary 
conditions show that the treatment of the boundary 
can have a significant effect on the first harmonic, 
unsteady pressure distributions for the cases with 
zero inter-blade phase angles. Using first or second- 
order temporal accuracy in the numerical algorithm 
did not make a significant difference in the unsteady 
pressure distributions for the present solutions. 
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FIGURE 1. - DEFORMING GRID TECHNIQUE, SIMPLIFIED GRID WITH 
EXAGGERATED MOTION. 
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FIGURE 1?. - MEAN FLOW PRESSURE DISTRIBUTION FOR BI- 
CONVEX AIRFOIL CASCADE, (5, = 60°, T S3 0 , g/c - 
0.767. T' - 0.07. 
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FIGURE 13. - FIRST HARMONIC PRESSURE DISTRIBUTION FOR BICONVEX AIRFOIL CASCADE, M 
pi = 60°+ 1 .2°, k = 0.185, 0 - -90°. 
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FIGURE 16. - BOUNDARY CONDITION STUDY FOR FIRST HARMONIC PRESSURE DISTRIBUTION FOR BICONVEX 
AIRFOIL CASCADE, M 1 = 0.80, = 60°±1.2°, k= 0.18S. 
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